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We have solved the quantum version of the Mattis model with infinite-range interactions. A 
variational approach gives the exact solution for the infinite-range system, in spite of the non- 
commutative nature of the quantum spin components; this implies that quantum effects are not 
predominant in determining the macroscopic properties of the system. Nevertheless, the model 
has a surprisingly rich phase behaviour, exhibiting phase diagrams with tricritical, three-phase 
and critical end points. 
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§1. Introduction 

Quantum spin systems with randomness are of ac- 
tive current interest because the interplay of randomness 
and quantum fluctuations often leads to nontrivial be- 
haviour. 1,2 J However, full exact analysis of such systems 
is very difficult, partly due to the randomness and partly 
due to non-commutativity of quantum spin operators. 
We therefore solve in the present paper a quantum ver- 
sion of the Mattis model with infinite-range interactions 
to investigate the effects of coexistence of randomness 
and quantum fluctuations. 

The Mattis model was originally proposed as a sim- 
ple spin glass model exhibiting only randomness, but no 
frustration. 3 ) It consists of Ising spins interacting via un- 
frustrated random exchange interactions. The random- 
ness can be gauged away, giving a simple ferromagnetic 
Ising model. The model is nevertheless nontrivial under 
an external field, in which case the problem reduces to 
that of the Ising ferromagnet with random local fields. 
In this paper we study the quantum version of the Mat- 
tis model in an external field. This is a more complex 
problem because, even in the absence of external fields, 
one cannot gauge away the randomness in the exchange 
interactions without violating the commutation relations 
of the quantum spin operators. 

There have been several investigations of the quantum 
Mattis model including the dispersion relation 4 ' and the 
symmetry of the ground-state. 5 ' These studies show that 
the non-commutativity of spin operators does not nec- 
essarily play an important role in the determination of 
the qualitative behaviour of the system. However, there 
has to date been no explicit solution for the equilibrium 
behaviour of the infinite-range model, even without an 
external field. The exact solution given in the present 
paper fills this gap; it also shows explicitly that, as was 
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to be expected from the earlier studies cited above, 4 ' 5 ' 
quantum fluctuations do not affect the structure of the 
phase diagram in an essential way. The physical reason 
for this fact is that the effective field acting on a spin 
is the sum of very many other spins when the range of 
interactions is infinite; but the sum of many quantum 
spin operators behaves like a classical vector. Hence, as 
we show, the problem reduces effectively to that of a 
single-site quantum spin in a classical external field. 

In spite of this fact, the system exhibits very rich phase 
behaviour. There are essentially three ordered phases, 
two with collinear spin orientations (parallel or antipar- 
allel to the external field), and one with non-collincar 
spins. These phases are separated by first or second or- 
der transition lines which terminate or meet at critical, 
tricritical, three-phase or critical end points. All these 
types of phase behaviour can be explicitly and exactly 
derived by relatively simple but nevertheless nontrivial 
manipulations. 

The paper is organized as follows. The model is de- 
fined and its variational free energy derived in §2. Ex- 
trcmization of the variational free energy gives the exact 
solution of the infinite-range model, even for quantum 
spins. (In an appendix, we also confirm the variational 
result by a direct explicit calculation.) The three special 
cases of the Ising, Heisenberg and XY models, which 
already exhibit the three types of ordered phases men- 
tioned above, are treated in §3. In §4, then, the model 
is analysed for general values of the anisotropy parame- 
ter in the exchange interaction. While no additional or- 
dered phases appear, the resulting phase diagrams now 
exhibit new and nontrivial features such as tricritical, 
three-phase and critical end points. The final section is 
devoted to conclusions. 

§2. Model and free energy 

The system we consider consists of N spin- 1/2 quan- 
tum spins Si. We use the rescaled spins erj = (2/%)Si, 
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whose components have eigenvalues <rf = ±1 (a = 
x,y,z), to describe the state of the system; this makes 
correspondences with classical Ising models easier to see. 
The Hamiltonian is defined as 



(Ti. 



(2.1) 



i<j 



Here the sum runs over all pairs of sites, J is a general 
3x3 coupling matrix, B is an external field, and the 
£j = ±1 are quenched random variables. If only J zz 
and B z are nonzero, then only the z-components of the 
spins, of, appear in H. Because the of all commute with 
each other, the quantum nature of the problem is then 
irrelevant, and one recovers the classical Mattis model 3 ' 
which is formulated in terms of Ising spins. 

To solve the model (2.1), it is tempting to try to gauge 
away the quenched disorder by the transformation er^ — > 
£i<7j. But this is impossible, because the gauged spins 
would no longer obey the required commutation relations 



(2.2) 



with e a( 3 7 the fully antisymmetric unit tensor. Instead, 
we solve the model using a variational mean field the- 
ory which treats the spins as uncorrelated with each 
other. Because we are dealing with a model with infi- 
nite range interactions, this approximation becomes ex- 
act in the thermodynamic limit TV — ► oo; an explicit self- 
consistency argument for this fact is given in Appendix 
A. A direct solution of a special case given in Appendix 
B also confirms the variational result. In the mean field 
approach, we start from a trial Hamiltonian 

H = — ^ hi ■ (Ti 

i 

with the associated variational free energy 

F = F + (H - H ) . 

Here (. . .) denotes an average over the Boltzmann distri- 
bution defined by H , and F is the corresponding free 
energy. The fact that mean field theory is exact then 
means that the minimum of F with respect to the hi 
is not just an upper bound on the true free energy F, 
but in fact equal to it in the thermodynamic limit. One 
easily evaluates F as 

F = -T^ln2cosh/3|/i i | 



N ^ 



^jirii ■ Jrrij - y^X B - h i) ■ m i> 



where (5 = 1/T as usual. The magnetizations are given 
by 

m l = (<7i> = t(/3hi) 

in terms of the "vector hyperbolic tangent" 

t(v) = ^-tanh(\v\). 
\v\ 



Minimizing F with respect to the hi gives the conditions 



X t 



h, B C-^-XX- 7 '" 



0, 



where the Xi are the local susceptibility tensors 

drrii 



Xi 



dh t 



(2.3) 



(2.4) 



with components given by xf 13 = dmf/dh^. It is easy 
to show that the Xi are positive definite and hence in- 
vertible; explicitly, one finds (hi = \hi\) 



hi <g> hi 



hi 



(1-tanh 2 (3hi)+( 1 - 



hi ® hi \ tanh j3hi 
phi 



with 1 denoting the 3x3 identity matrix. Hence the 
minimization conditions (2.3) can be written as 

hi = B + ti^J2tj Jm 3 



= B + £iJ(n + m + — n_m_) — — Jrrii, (2.5) 



where 



m± 



Nn-i 



E 

iei± 



are the magnetizations of the sublattices I± = {i : ^ = 
±1}, respectively; n + and n_ give the fraction of spins 
contained in the two sublattices. The last contribution to 
hi in (2.5) can be neglected in the thermodynamic limit, 
and so the local variational fields become dependent only 
on the sublatticc, not the site i itself: all the spins in 
the same sublattice feel the same local field and hence 
have the same magnetization. One can now rewrite the 
variational free energy per spin, / = F/N, in terms of 
the sublattice magnetizations alone, with the result 

/(m+,ra_) = — Tn + s(\m + \) — Tn_s(|m_|) 

— i(n + m + — n_m_) • J(n + m + — n_m_) 

— B ■ (n + m + + n_m_), (2.6) 

where the entropic contribution is expressed as usual in 
terms of the entropy of a binary distribution, 



s(m) = — 



1 



m 1 
■ m ■ 



m 



1 



m , 1 
■ In ■ 



m 



(2.7) 



2 2 2 2 

The true free energy per spin is obtained as the minimum 
of the variational free energy, hence the final result 



/ = m^_ /(m+ ' m - ) - 



(2.8) 



In the following, because the variational free energy is ex- 
act for N — > oo, we drop the tilde on /. The requirement 
that f(m + , m_) must be stationary with respect to m + 
and ra_ gives two self-consistency equations which can 
be written in the compact (and intuitively appealing) 
form 

m + =t[(n + Jm + -n-Jm- + B)/T] (2.9) 
m-=t[(-n + Jm + +n-Jm-+B)/T] (2.10) 
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We emphasize that the energetic terms in the free en- 
ergy (2.6) are exactly of the form that one would expect 
if the sublattice spins 

<T± = T7— E (Tl 



Nn± 



iei± 



were classical vectors rather than quantum operators. 
Working out commutation relations such as 



+ iel+ 7 



2i 



E e "/37 cr + 



indeed shows that in the thermodynamic limit N — > oo 
the vectors cr± become classical, with all components 
commuting with each other. Note that even though the 
m± are classical, the quantum (spin) nature of the prob- 
lem is still reflected in the functional form (2.7) of the 
cntropic contribution to the free energy (2.6). 

Having derived the free energy of our model for general 
J and B, we specialize from now on to the case 



J = diag(l,l,A), B = (0,0, B). 



(2.11) 



The parameter A here interpolates between three im- 
portant limits: For A = and 1, respectively, we have 
the XY and Heisenberg versions of the model, while for 
A — > oo the classical (Ising) Mattis model is recovered. 
We will analyse these three cases separately first before 
studying the richer behaviour obtained for intermediate 
values of A. Because all phase diagrams are symmetric 
under B — > — B, we generally restrict ourselves to B > 0. 

Finally, before proceeding, we note that with the 
choice (2.11), the model has a trivial rotational sym- 
metry in the xy-planc. We break this symmetry by re- 
quiring that m!i_ = and m x + > 0. Minimization of (2.6) 
with respect to rotations of m_ implies that the com- 
ponents of m + and m in the a;y-plane are antiparallel 
to each other, thus also m v _ = 0, and m x _ < 0. So the 
minimization in (2.8) only has to be carried out over 
the four magnetization components m z +1 mi, m x _, m x _, 
with the latter being respectively non-negative and non- 
positive. We also assume throughout that the positive 
sublattice J + contains more spins than the negative one, 
i.e., n + > n_, and take n_ > in order to exclude 
the trivial case of an anisotropic quantum ferromagnet 
without disorder. Instead of n + and n_, we will some- 
times use the parameter e, defined by n± = (1 ± e)/2; 
e — ► 1 then corresponds to the disorder- free limit, and 
our assumptions on n+ and n_ translate into < e < 1. 

§3. Ising, Heisenberg and XY models 

3.1 Ising limit (A — > oo^ 

In the limit A — > oo, the free energy (2.6) is min- 
imized when the sublattice magnetizations point along 
the z-axis. We thus only have two nonzero order pa- 
rameters m z + and ml to consider, and effectively re- 



cover the Mattis model with classical Ising spins. Using 
s(m) = s(—m), the free energy simplifies to 

//A = -fn+s(m z + ) - fns(m z _) 

- ^(n + m z + - n_m z _) 2 - B(n + m z + + n_m z _), (3.1) 

and the stationarity conditions (2.9) and (2.10) become 

m* = tanh[(n+m^ - n^ra z _ + B)/f] (3.2) 

m z _=t&nh[(-n + m z + +n-m z _+B)/f]. (3.3) 

Here we have introduced the rescaled temperature T = 
T/A and field B = T/A which are the relevant control 
parameters for A — > oo. 

Considering first the zero temperature limit T — > 0, 
the only possible solutions of (3.2) and (3.3) are m z ± — 
±1, so we only need to compare the four resulting val- 
ues of /. For large B, one finds in this way that 
m z + = m z _ = 1. We call this phase A + , to indicate that 
both sublattice magnetizations m + and ra_ are aligned 
along the direction of the field B, with both pointing in 
the same direction. At B = n+, there is a first order 
transition to m z + = 1, m z _ = —1; we denote this new 
phase A_ because the sublattice magnetizations, while 
still aligned with the field, now point in opposite direc- 
tions. Finally, at B — we have the conventional first 
order transition where both magnetizations change sign, 
and then the mirror image of the A + -A_ transition oc- 
curring at B = — n + . As pointed out earlier, all phase di- 
agrams of our model have this symmetry under B — > —B, 
so we will not mention results for B < in the following. 

In the T-B phase diagram, the first order transitions 
found above for T = mark the beginnings of two first 
order transition lines, both ending in critical points. The 
general procedure by which we find such points is out- 
lined in Appendix D; after a little algebra, the two rel- 
evant conditions (D-2) and (D-5) become 

f = n+{l - (m z + ) 2 } + n_{l - (m z _) 2 } 

= n+m z + [l - (m z + ) 2 } - n_ml[l - (m z _) 2 }. (3.4) 

These need to be solved along with (3.2) and (3.3). The 
critical point that marks the end of the first order transi- 
tion line at B = has m z + = m z _ = and is thus located 
at T = n + + ra_ = 1, B = 0. The critical point ter- 
minating the A_|_-A_ transition line, on the other hand, 
has to be found numerically, along with the location of 
that line itself. The resulting phase diagram is shown in 
Fig. 1 for e = 0.01. 

Note that because of the critical point terminating 
the line of first order transitions between A + and A_, 
these two phases can be smoothly transformed into one 
another by moving along a continuous path in the T- 
B phase diagram. Along this path, the value of m z _ 
changes sign by passing through zero; the points where 
this happens obey B = n + tanh(2_B/T), as is easily de- 
rived from (3.2) and (3.3). As in the case of the liquid- 
gas transition, this makes the thermodynamic distinction 
between the two phases somewhat arbitrary. Unlike the 
traditional case, however, the critical point is not due to 
a symmetry under sign reversal of the magnetization; the 
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Fig. 1. Phase diagram for the Ising limit A — > oo, for e = 0.01. 
The dashed lines indicate first order transitions, and end in crit- 
ical points. The transition at B/A = and T/A < 1 is a 
conventional one, where all magnetizations change sign; at the 
other transition (B/A ^ 0), the relative orientation of the sub- 
lattice magnetizations changes from antiparallcl (A_) to parallel 
(A_|_). The dotted line indicates the points where the A_ and 
A+ phases transform smoothly into each other as mf_ passes 
through zero. As the inset shows, the A_— A+ critical point is 
not on this line; it has m z _ > 0. 



critical value of m z _ is therefore nonzero (and, by (3.4), 
positive). This implies that just below the critical tem- 
perature, the A_|_-A_ transition is actually between two 
phases with mi > 0, rather than two phases with oppo- 
site signs of m z _ . 

The case of e = in the Ising limit is worth an atten- 
tion although we generally assume e > 0. When e = 0, 
the system is equivalent to the random-field Ising model 
with symmetric distribution of randomness. This prob- 
lem has already been solved by the mean-field approxi- 
mation. 6 ) The phase diagram is similar to Fig. 1 except 
that the smooth crossover between the A + and A_ phases 
along the dotted line is now replaced by a second order 
transition. This fact can be verified easily, for example, 
in the case of B = by setting n + = ra_ = 1/2 in (3.2) 
and (3.3) and rewriting these equations for the combina- 
tion m z + — m z _ . 

Finally, we note that the results (3.1,3.2,3.3) are valid 
not only in the limit A — > oo, but for all A > 0, as long as 
both sublattice magnetizations m + and m_ point along 
the z-axis. This means that the properties of the A + 
and A_ phases are independent of A when expressed 
as functions of B and T. In particular, the first order 
A + A_ transition as well as the first order transition at 
B = 0, T < 1 will be present in all phase diagrams unless 
"masked" by other phases. 

3.2 Heisenberg model (A = l) 

In the Heisenberg case A = 1, the coupling ma- 
trix (2.11) is isotropic, and the free energy (2.6) can be 
rewritten as 

/ = — Tn + s(m + ) — Tn_s(m_) + ^to 2 

-n^m^-ii 2 m 2 _-5m 2 , (3.5) 



where m± = \m±\, and to and m z denote the modulus 
and z-component, respectively, of the average magneti- 
zation 



rn 



n_ra_ 



of all spins. The only dependence on the orientation 
of the magnetizations is through the last term, which 
takes its minimal value —Bm (remember that we assume 
B > 0) when m z = to, i.e., when m is parallel to B. 
This gives the free energy 



/ = —Tn + s(m + ) — Tn. 



2 2 

-n , m , 



2 2 

n_m_ 



_s(m_) + -^m 2 
Bm 



(3.6) 



which needs to be minimized subject to the constraints 
\n + m + — n_m_| < m < n + m + + n-in-. Note that the 
values of m + , m_ and to determine the orientations of 
the magnetizations uniquely. This is clear geometrically 
because n + m + , n^m^ and m all lie in the same plane 



(the xz-plane, by our assumption that rrv\ 



0) 



and form a rigid parallelogram; see Fig. 2. The angle of 
a possible rotation of this parallelogram about the origin 
is fixed by the requirement that m must be parallel to B 
(i.e., along the positive z-axis for B > 0). The remaining 
indeterminacy with respect to a reflection about the z- 
a xis is removed by our assumption that m x + > 0. 



n m 




n in 




Fig. 2. Geometry of the magnetizations for the Heisenberg 
model. Shown are (in the xy-plane) the total magnetization m = 
n+Tn.+ + n_m_ and the sublattice magnetizations scaled by 
the fractions of spins in the two sublatticcs, ji+tti+ and n— tyi— . 
The moduli of the sublattice magnetizations m+ = \m+\ and 
m_ = |m_| are kept constant while m = \m\ is decreased from 
left to right. For the maximum value of m, we have an A_|_ 
phasc (left), then we obtain an R phase, with m_ rotating from 
the positive (upward) towards the negative (downward) z-axis. 
Finally, when m assumes its minimal value (right), we obtain 
an A_ phase. This kind of sequence is observed when pass- 
ing through the R phase in the Heisenberg phase diagram by 
decreasing the field B at constant T; the moduli of the sublat- 
tice magnetizations remain constant, and they rotate in such a 
way as to keep m pointing along the direction of the field (the 
z-dircction). 



Focussing now on the minimization of (3.6) with re- 
spect to m, we note that for B > n^m + + n_m_ this 
minimum occurs at the maximum value of to, to = 
n + m + + n_TO_. Geometrically, this implies that both 
to + and to_ are directed along the positive z-axis, so we 
have an A + phase. The free energy becomes 

/ = —Tn + s(m + ) - Tn_s(m_) - i(n + TO + - n_m_) 2 



-B(n + m + + n_m_), 
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and the sublattice magnetizations therefore obey 

m+ = tanh[(n + m + - ra_m_ + B)/T] (3.7) 

m_ = tanh[(-n+m+ + n_m_ + B)/T]. (3.8) 

For small fields obeying B < n + m + — n_m_, on the 
other hand, / is minimized at the minimum value of m, 
m, = n + m + — n_m_. With m now pointing along the 
negative z-axis, we have an A_ phase with free energy 

/ = — Tn + s(m + ) — Tn_s(m_) — -(n+ra+ + n_m_) 2 

—B(n+m + — n_m_), 

and corresponding self-consistency equations 

m + = tanh[(n + m + + n_m_ + B)/T] (3.9) 

m_ = tanh[(n+m+ +n_m_ - B)/T]. (3.10) 

The above results are identical to (3.1), (3.2) and (3.3) for 
the Ising case, bearing in mind that A = 1 and 777+ = m z + 
here, and that m z _ = ±m_ in the A + and A_ phases, re- 
spectively. This conclusion agrees with our general state- 
ment in § 3.1 that properties of the A + and A_ phases 
are independent of A. 

For intermediate values of B, finally, the minimum of 
the free energy (3.6) occurs at a non-extremal value of 
m. In this regime, we have a new, rotated (R) phase 
where neither m+ nor m point along the z-axis. Min- 
imization of (3.6) with respect to m gives now m = B, 
and thus 

/ = — Tn + s(m + ) — Tn_s(m_) — n 2 + m 2 + — n 2 _m 2 _ — ^B 2 . 

Stationarity with respect to m+ and to_ yields the self- 
consistency conditions 



m + = tanh(2n+m+/T) 
m_ = tanh(2n_m_/T) 



(3.11) 



in which the moduli of the sublattice magnetizations are 
decoupled and the field B no longer appears. This may 
seem surprising at first, but has a simple explanation: 
From (2.9) and (2.10), the effective fields that determine 
the sublattice magnetizations are (using J = 1 for A = 
1) 



h± = ±(n+ra+ — n_m_ 
But in the R phase, m = _B, so that 



B. 



h± = ±(n+ra+ — n_m_) + (n + m + + n_m_) = 2n±m± 

and the coupling of the sublattices and dependence on 
B disappear, in agreement with (3.11). The geometrical 
implication of (3.11) is that only the orientations but not 
the moduli of the sublattice magnetizations change as 
the R phase is traversed. Figure 2 shows that as m (and 
thus B) decreases, m_ rotates away from the positive 
z-axis and towards the negative z-axis; when it reaches 
the latter, a transition to an A_ phase occurs. The mag- 
netization of the positive sublattice, m+, first rotates 
away from the z-axis and then back towards it; this fol- 
lows from the fact that m = n + m + + n^ni- must keep 
pointing along the positive z-axis. 



Above, the A+-R and R-A_ phase boundaries were 
given as B = n+m+ ± n_m_, respectively. At T = 0, 
where m + = to_ = 1, these reduce to B = n + + n_ = 1 
and B = n + — ri- = e; for nonzero temperatures, we have 
to find the phase boundaries numerically by solving the 
relevant self-consistency equations for m, + and m_. We 
then find (see Fig. 3) that there is a line of second order 
transitions connecting the zero-temperature A + -R and 
R-A_ transitions; the R phase only occurs inside the 
"loop" formed by this line. The only other feature of the 
phase diagram is the first-order transition line at B = 0, 
T < 1, which is identical to the one found for the Ising 
case. 

1.2 i ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 




0.2 0.4 



0.6 
T 



0.8 



1.2 



Fig. 3. Phase diagram for the Heisenberg case A = 1, for e = 0.4. 
The dashed line shows the conventional first order transition at 
B = 0, where all magnetizations change sign; the solid line is a 
line of second order transitions between aligned (A_|_ or A_) and 
rotated (R) phases. The sublattice magnetization m_ passes 
smoothly through zero on the dotted line. 



3.3 XY model (A = 0) 

As the final "simple" case we consider the XY ver- 
sion of our model, which is obtained for A = 0. In 
this case it is convenient to introduce the mirror im- 
age (in the xy-plane) of m about the z-axis, given 
by m_ = (— m x _, 0, mi), and the corresponding total 
(pseudo-) magnetization rh = n + m + + n^rh-. In these 
variables, the free energy (2.6) becomes 

/ = — Tn + s(m + ) — Tn_s(rn_) — — Bfh z 



= —Tn + s(m + ) — Tn_s(rn_) 

1 2 1 2 

--m +-m z -Bm z . 



(3.12) 



Minimizing over the allowed values of m z £ [— m, m] 
gives tti z — B as long as B < rh. The last two 
terms in (3.12) then reduce to the constant —B 2 /2, 
and the remainder of the free energy is minimized (for 
given TO+ and rh-) when 777 takes its maximum value 
777 = 77+777+ + 77_7n_. Geometrically, this means that 
m+, rh and rh are all parallel to each other. The free 
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energy is 

/ = -Tn + s(m + ) - Tn_s(ra_) 
~(n + m + + n_rh_) 2 - ^B 2 , 
while the stationarity conditions 

m + = tanh[(n + m + + ra_m_)/T] 
m_ = tanh[(n + m + + ra_m_)/T] 

show that m + = m_, hence also rh = m + = rh- 
with rh = tanh(ra/T). The three vectors m + , rh and 
rh are therefore not just parallel, but in fact identical; 
their orientation in the xz-plane is given by the ratio of 
rh z = B and rh. Reverting to the original vectors, we 
have that m + and m_ are rotated away from the z-axis 
and are mirror images of each other under a reflection 
about this axis; we therefore have an R phase. The av- 
erage magnetization m always points along the positive 
z-axis. Starting from B = 0, ra_ and m + are directed 
along the negative and positive ai-axes, respectively (this 
follows from rh z = and m5_ > 0). As B is increased, 
both sublattice magnetizations rotate towards the z-axis 
which they reach at the point where B = rh. 

For larger B, we have an A+ phase. In this phase, 
rh z = rh; inserting this into (3.12), the minimum with 
respect to rh is again reached for rh = n + m + + n_m_, 
giving 

/ = — Tn + s(m + ) — Tns{fh-) — B(n + m + + n_m_), 

and thus m + = rh- = rh = tanh(B/T). All vectors m + , 
rh and rh are again identical to each other; because 
they are now oriented along the z-axis, the same is true of 
the original magnetizations m + , m_ and m. In Fig. 4, 
we show the resulting phase diagram; as expected, there 
is a line of second order transitions between the R and 
A+ phases. Note that obtaining this line numerically is 
trivial for the XY model: Combining rh = tanh(m/T) 
and B = rh gives B = tanh(S/T). It follows in particu- 
lar that in the XY limit, the phase diagram is actually 
independent of the amount of disorder (as specified by 

It may appear strange that the disorder has no effect 
on phase behaviour in the XY limit A = 0. But there is 
in fact again a very simple explanation for this. In the 
model (2.1), make a gauge transformation only on the x 
and y-components of the spins, of — > £jcf, erf — > £jof . 
This leaves the commutation relations (2.2) unaffected; 
indeed, it is just a rotation by tt around the z-axis of 
the spins with £ = — 1. But this transformation actu- 
ally gauges away the disorder completely, so all results 
must be independent of e, in agreement with our findings 
above. Note that this argument relies on the fact that 
only the x- and y-components of the spins appear in the 
disordered (interaction) part of the Hamiltonian (2.1), 
and is therefore restricted to A = 0. 

§4. Phase behaviour for general A 

We now turn to the study of our model (2.1) for gen- 
eral values of A. The differences in the phase diagrams 
for the three cases studied above (A = 0, 1 and oo) 




0.2 0.4 0.6 0.8 1 1.2 
T 



Fig. 4. Phase diagram for the XY case A = 0. The solid line 
indicates the second order transition between the rotated (Ft) 
and the aligned (A_)_) phases. Note that this phase diagram is 
independent of e, i.e., of the amount of disorder. 



already suggest that nontrivial phase behaviour may oc- 
cur for intermediate values of A. To orient ourselves, we 
consider first the zero temperature limit, using A and B 
as the axes of our phase diagram (and considering e as 
fixed) . 

For large A, we expect essentially Ising behaviour, 
with an A + phase for large B, a first order transition to 
A_ at B = n + A, and the conventional first order transi- 
tion at B = where all magnetizations change sign. In 
fact, as explained in §3.1, this will be the zero tempera- 
ture phase behaviour for general A unless other phases 
intervene. We therefore study next the limits of stabil- 
ity (i.e., the spinodals) of the A + and A_ phases with 
respect to a transition to the R phase. This is easiest if 
we write the sublattice magnetizations as 

m + = m + (sin + , 0, cos <j) + ) , 

m = rri- (— sin <j>_, 0, cos </>_). 

Here 4>± are the angles (in the rrz-plane) that m + and 
m make with the z-axis; the signs of the angles were 
chosen such that our conventions m 3 ^ > and m x _ < 
always imply non-negative angles. Using that m + = 
rri- = 1 at T = 0, the free energy (2.6) thus simplifies 
to 

/ = — i(n + sin0 + + ra_ sin</>_) 2 

— — A(n + cos — n_ cos (f>-) 2 

—B(n + cos (f> + + n- cos (4.1) 

With only two order parameters (<p + and 0_) remaining, 
it is straightforward to find the matrix of second deriva- 
tives of /. The criterion (D-2) implies that a spinodal 
instability occurs when the determinant of this matrix 
vanishes. Evaluating the latter for <p + = <p- =0, one 
finds for the spinodal of the A + phase the condition 

B 2 — B + e 2 (A — A 2 ) = 0, (4.2) 
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while for the A_ phase (<p + = 0, <fi- = n) the corre- 
sponding result is 

- J B 2 + J Be + A 2 -A = 0. (4.3) 

The signs of the expressions on the l.h.s. have been 
chosen such that they are positive when the phases 
are stable. Also, because of the symmetry of the free 
energy (4.1) under <fi + — > — <p + and 0_ — > — or 
— > 2ir — 0_ , these spinodals automatically satisfy the 
critical point criterion (D-5) and so are in fact critical 
points. Figure 5 shows a plot of the spinodals lines (4.2) 
and (4.3) for e = 0.4. We see that the A + phase is stable 
for large fields, but destabilizes as B is lowered; the A_ 
phase, on the other hand, tends to be stable for smaller 
values of B, and large A. Nontrivially, however, the A_ 
spinodal shows re-entrance: For A € [A rc , 1], with 

A rc = i(l + V / l-e 2 ), 

Eq. (4.3) has two physical solutions for B, so the A_ 
phase is unstable at zero field, becomes stable at in- 
termediate values of B, and destabilizes again as B is 
increased further. 



CQ 1 
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re/ tn 
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A_ 

Ising 



0.5 1 1.5 2 2.5 
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Fig. 5. Zero temperature phase diagram for e = 0.4. The transi- 
tions between the three phases A+, A- and R arc indicated by 
bold lines (solid for second order, dashed for first order). Dot- 
ted lines show the continuations of the phase boundaries into 
mctastablc or unstable regimes. Note that the A_-R transition 
is re-entrant for A ro < A < 1. Also, for A tr i < A < Ai s j ng , the 
first order A_|_-A_ transition is — because of the instability of 
the A_|_ phase — pre-empted by a first order transition from A_ 
to R. Below A tr ;, this transition is second order, implying that 
at A = A tr i there is a tricritical point (marked by the circle). 



In Fig. 5, the first order A + -A_ transition line B = 
n + A is also shown. Moving along this transition line 
from large to small A, the first spinodal which one crosses 
is that of the A + phase, at the value of A given by 

A -g 2e+1 
Ising 3e + 1 " 

For large values of A, the instabilities of the A + and 
A_ phases are therefore pre-empted by the first order 
A_|_-A_ transition. In this regime, we expect pure Ising 
behaviour even at non-zero temperature, and this is in- 



deed what we find (see below). 

Now consider a value of A just below Ai S i ng . The 
A_ phase is stable for B = 0, and remains so until 
at B = n + A the free energy of the A + phase becomes 
lower. The latter is still unstable, however, because we 
are below the A + spinodal. There must therefore be 
a stable phase with lower free energy. This phase can 
only be an R phase (it is neither A + nor A_), and a 
first order transition to this phase must actually occur 
at B < n + A. This implies that there is a line of first or- 
der A_-R transitions extending to the left of the point 
A = Aisi ng ,_B = n+Aising. Where this line meets the 
A_ spinodal (i.e., the line of second order transitions 
between A__ and R), there will be a tricritical point. Ap- 
plying the criterion (D-6) to the free energy (4.1), one 
finds after some algebra that this point obeys, in addi- 
tion to (4.3), 

B 2 = A(A-l)(4A-3). 

It can be shown that, as e varies between and 1, this 
tricritical point moves smoothly from A = 1 , B = to 
A = 3/2, B = 3/2. In particular, if we call A tr ; the value 
of A at the tricritical point, we have 1 < A tr ; < 3/2 < 
Ai s i n g < 2 for all e. 

Having clarified the structure of the zero temperature 
phase diagram, we can now move on to the finite tem- 
perature case. The numerical results we show were ob- 
tained as follows: For spinodals and tricritical points, 
we derived analytically the form of the relevant condi- 
tions (D-2) and (D-6) (for our free energy (2.6) with the 
four order parameters m z + , m z _, m x + and m x _). We then 
solved these numerically, along with the self-consistency 
equations (2.9). First order transitions were located 
as usual by comparing the free energies of the relevant 
phases. All results were obtained from double-precision 
routines for solving nonlinear simultaneous equations, 
and cross-checked using a symbolic manipulation soft- 
ware package with "arbitrary precision" floating point 
operations. We can distinguish a total of seven different 
phase diagram topologies, depending on the value of A: 

Regime 1: XY -like behaviour (0 < A < A ro ). For 
small values of A, the phase diagram has essentially the 
same features as for the XY limit A — > 0; an example 
is shown in Fig. 6. It is clear that this behaviour cannot 
persist up to A = 1, however: From the zero tempera- 
ture phase diagram, we know that for A ro < A < 1 , there 
must be re-entrant behaviour. Correspondingly, the sec- 
ond order R-A + transition line must develop a "dent" — 
as if someone was pushing against it from the positive 
T-direction — as A increases; this dent will reach T = 
exactly at A = A rc . Before this happens, re-entrance will 
already be visible for nonzero T; Fig. 6 confirms this. In 
principle, one could use the appearance of re-entrance 
at nonzero T to further divide this regime into two sub- 
regimes, but we choose not to do so because the overall 
topology of the phase diagram remains unchanged. 

Regime 2: Between XY and Heisenberg (A ro < A < 
1). At A = A ro , the R-A + transition line "pinches off" 
at the -B-axis; for larger A, we therefore have two sep- 
arate transition lines of this kind (see Fig. 7). For low 
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cq 




Fig. 7. Phase diagram in regime 2 (A ro < A < 1), for e = 0.4 
(where A rc = 0.958 . . .) and A = 0.96. The lines have the same 
meaning as in Fig. 6. Re-entrant behaviour now occurs even at 
zero temperature. 



(see Fig. 8). At A = A tr i, the tricritical point that we 
found earlier appears on the T — axis, marking the 
transition to the next regime. 
1 r- 



1 ■ 1 ■ ^ — 1 1 

0.5 1 1.5 

T 

Fig. 6. Phase diagrams in regime 1 (A < A ro ) for e = 0.4. The 
solid lines mark second order transitions between the phases. As 
in Fig. 1, the dotted lines show where the A_ and A+ phases 
transform smoothly into each other (mf_ = 0). For small A (top, 
A = 0.9), the phase diagram resembles qualitatively that of the 
XY model (Fig. 4). In the bottom graph, A = 0.95 is close to 
A rc = 0.958..., and re-entrant behaviour appears at nonzero 
temperature. This is a precursor of the transition to regime 2. 



enough temperatures, the sequence of phases observed as 
B is increased from is therefore R-A_-R-A + . As the 
Heisenberg case A = 1 is approached, the R-A_ transi- 
tion line (a loop, if we bear in mind its mirror image for 
negative B) moves closer to the horizontal axis B = 
and at A = 1 collapses into the first order transition line 
of Fig. 3. The latter exists for all A > 1 (extending up 
to T = 1, i.e., T = A, as expected from our discussion 
of the Ising case) and will not be mentioned explicitly in 
the following. 

Regime 3: Heisenberg-like (1 < A < A tri ). To make 
the eventual connection with the Ising limit A — > oo 
more apparent, we will use the rescaled variables T — 
T/A and B = B/A as the axes of all phase diagrams 
from now on. In the numerical examples, we also switch 
from e = 0.4 to e = 0.01, where the intermediate regimes 
explained below are somewhat easier to visualize. In the 
regime 1 < A < A tri , the phase diagram has qualita- 
tively the same shape as for the Heisenberg case A = 1: 
A loop of second order transitions between R and A + 
or A_, respectively, beginning and ending on the B-axis 




o 1 — — ----- — 

0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
T/A 



Fig. 8. Phase diagrams in regimes 3—7 (A = 1, 1.2, 1.4, 1.5, 2 
from the outside to the inside), for e = 0.01. Solid and dashed 
lines mark second and first order transitions, respectively; the 
circle indicates a tricritical point. See following figures for de- 
tails. 

Regime 4: Tricritical (A tr ; < A < A3). As A in- 
creases, the tricritical point moves out to larger T. To 
the left of it, the A_-R transition is now first order (see 
Fig. 8). One might naively expect that as A is increased 
further towards the Ising limit of large A, the "loop" 
of transitions between R and A_ or A + would simply 
collapse at A ls i ng onto the Ising first order A + -A_ tran- 
sition line. In fact, two other regimes appear first. 

Regime 5: Tricritical plus three-phase (A3 < A < 
Aqep)- It turns out that the R-A transition loop with 
the tricritical point on it does not just shrink along the B 
direction, but also moves towards smaller values of T as 
A is increased. Eventually, at some A3, the Ising A + -A_ 
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transition line therefore "pokes" through the loop. This 
must happen at a point where the transition is first or- 
der; at a point with a second order transition this would 
be impossible, as it would imply a first order transition 
between two phases which are actually identical, being 
both in critical coexistence with the same third phase. 
This means that one has a point where three phase tran- 
sition lines meet, i.e., a three-phase point; see Fig. 9. 
There is still a tricritical point where the R-A + phase 
transition changes from second to first order. Note that 
this regime generally corresponds to only a very small 
range of A; for e = 0.01, we estimate A3 w 1.39 and 
A CE p » 1.41. 
0.467 



cq 



0.466 
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0.464 
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A_ 





0.584 0.586 0.588 0.59 
T/A 



0.592 



Fig. 9. Detail of a phase diagram in regime 5 (A = 1.4, e = 0.01). 
Solid and dashed lines mark second and first order transitions, 
respectively; circles indicate tricritical points. There is a three- 
phase point where the three first order lines meet. 



Regime 6: Critical end point (Acep < A < Ai s i ng ). 
The relative positions of the tricritical and three-phase 
points in the previous regime change as A increases, un- 
til at some Acep they coincide. From there onwards, 
the tricritical point is no longer accessible (it is in a 
metastable or unstable part of the phase diagram). In- 
stead, as shown in Fig. 10, one now has a line of second 
order A + -R transitions that meets a line of first order 
transitions (between A_ and R for small T, and between 
A_ and A + for large T) at a critical end point. As A in- 
creases, this point shifts towards lower temperatures and 
eventually disappears at A — Ai s ; ng . Fig. 11 illustrates 
the transition between regimes 5 and 6 by showing the 
curve traversed by the tricritical point as A is varied, to- 
gether with the Ising A + -A_ transition line. The value 
of A where these two curves meet (i.e., where the tri- 
critical point "collides" with the A + -A_ transition and 
thus turns into a critical end point) defines Acep- 

Regime 7: Ising-like (A > Ai s ; ng ). Finally, for large 
values of A, one has pure Ising behaviour, with a A- 
independent phase diagram (when represented in terms 
of B and T) exhibiting the by now familiar line of first 
order A + -A_ transitions (see Fig. 1). 
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T/A 



Fig. 10. Detail of a phase diagram in regime 6 (A = 1.5, e = 
0.01). Solid and dashed lines mark second and first order tran- 
sitions, respectively. There is a critical end point where the line 
of second order transitions meets the first order line. 
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0.464 



0.462 



0.62 



Fig. 11. The transition between regimes 5 and 6, for e = 0.01. 
The dashed line shows the Ising A+-A_ transition line (which, 
when plotted in terms of B = B/A and T = T/A as done here, 
is independent of A). The solid line shows the curve traversed 
by the tricritical point as A is varied, moving from right to left 
with increasing A. The crossing of the two curves (i.e., the value 
of A for which the tricritical point meets the A+— A— transition 
and thus turns into a critical end point) defines Acep- Note 
that for lower values of A than shown here, the tricritical point 
would first continue to move right, but eventually swing back 
towards the S-axis, meeting the latter at A = A tr i- 



§5. Conclusion 

We have solved the infinite-range quantum Mattis 
model by a variational method that gives the exact so- 
lution in the thermodynamic limit. The model has var- 
ious interesting aspects such as randomness (although 
without frustration), quantum effects and competition 
between exchange interaction and external field. We 
found three ordered phases, two of which have spin states 
collinear with the external field and the remaining one 
with non-collinear rotated spin states. The phase dia- 
gram has a very rich structure depending upon the var- 
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ious parameters, in particular the anisotropy of the in- 
teraction. 

We now ask how important quantum effects are in pro- 
ducing the intricate macroscopic behaviour of the sys- 
tem. The form of the entropy term in the free energy 
(2.6) is a direct consequence of the spin-1/2 character- 
istics of a single quantum spin. The energy term, on 
the other hand, is of the form one would expect classi- 
cally. Thus it is clear that the T = properties do not 
reflect quantum effects. Finite-temperature behaviour, 
on the other hand, could be affected by quantum fluc- 
tuations. This fact is somewhat counter-intuitive since 
quantum effects are usually most important at low rather 
than high temperatures. An interesting manifestation of 
finite-temperature quantum effects is seen in the equa- 
tions of state (3.8) to (3.10) of the Heisenberg model in 
the A + and A_ phases. These equations of state have 
exactly the same form as for the Ising case (3.2) and 
(3.3), which would not be the case if the spins in the 
Heisenberg model were classical vectors. The transverse 
[i.e. x and y) components of the spin-1/2 operator have 
disappeared and only the z component comes into play, 
as would be expected for a single quantum spin in a clas- 
sical field. 

We stress that, even though the quantum nature of 
the individual spins in the model does manifest itself at 
finite temperatures, this does not produce any qualita- 
tively new ordered phases beyond those already found 
at T = 0. The possible types of ordering are thus de- 
termined by the classical considerations for the ground 
state; quantum effects only control how these ordered 
phases and the boundaries between them arc arranged 
in the finite temperature phase diagrams. 

Finally, a natural question that arises is how our re- 
sults would change in the more physically realistic case of 
finite range interactions (corresponding to finite dimen- 
sionality of the system). While we cannot give a definite 
answer to this question at this point, existing investiga- 
tions of finite-dimensional quantum Mattis model 4, 5 ) do 
suggest that the classical picture gives reasonable pre- 
dictions for some features of the model. On the basis of 
these results, we conjecture that our observations for the 
infinite-range model should give an at least qualitatively 
reliable guide to the finite-dimensional problem. Future 
work is obviously needed to clarify this point. 
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Appendix A: Self-consistency of mean field ap- 
proximation 

In this appendix, we show that mean field theory gives 
exact results for the generalized Mattis model (2.1). The 
intuitive reason for this is clear: The field that each spin 
experiences is an average of O(N) other spins and there- 
fore becomes nonfiuctuating (and classical) in the ther- 
modynamic limit. 



To see more explicitly why the mean field approxima- 
tion, which assumes that all spins are uncorrelated with 
each other, becomes self-consistent for N — > oo, we con- 
sider a generalization of our original Hamiltonian (2.1) 
to site-dependent fields: 

h = -— y^&gjo-j ■ j(Tj - y~] Bj ■ <ji. 

i<j i 

The fluctuation-dissipation theorem then relates the sus- 
ceptibility matrix to the spin-spin correlations according 
to 

d (<Tj) dF / \ \ \ \ , > 

M = -dB~dB- = * < ACT ' A ^ ■ = - • 

(AT) 

For given i and j, this is an equality between 3x3 tensors, 
whose components are written explicitly as d (erf) /dB 1 - 
= (3 (Acrf AcrJ). We can thus verify that mean field the- 
ory is self-consistent by working out the susceptibility 
matrix and using it to show that correlations between 
different spins vanish for N — > oo. 

To obtain the susceptibility matrix, we start from the 
mean field equations. By a direct generalization of (2.5), 
these are 

hi = B l +£, i J(n + m + -n-m-)- -^Jm,i. (A-2) 

For the case of site- independent fields Bi = B, they re- 
duce to (2.5), with the fields hi and magnetizations m,i 
being the same for all spins in each of the two sublat- 
tices I± (in the thermodynamic limit). We now add a 
small perturbation SB to the field of one of the spins; 
without loss of generality, this spin can be taken as <T\. 
For definiteness, we also assume that <J\ is in the posi- 
tive sublattice J + ; the calculation for the opposite case 
is completely analogous. The solution to the mean field 
equations (A-2) will then be such that all spins in the 
sublattice 7_ still have the same fields h and magne- 
tizations m . In I + , on the other hand, we have to 
distinguish between the field and magnetization of the 
chosen spin (hi and mi) and those of all other spins in 
the sublattice; we denote the latter by and m-i . Using 
that the average magnetization of this sublattice is now 

m+ = -J— [mi + (Nn + - l)m 2 ] 
the mean field equations (A-2) then take the form 



Subtracting the corresponding equations for SB = 
gives relations between the deviations of the fields (Shi 
etc) and magnetizations (Smi etc) from their values for 
site-independent fields: 



For small SB, we can linearize these using 
Smi = x+Shi, 5m,2 = x + <5h.2, 5m- = x~Sh-, 
where x+ an d X- are the local susceptibility ten- 
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hi = B + SB + J 
h 2 = B + J 
h =B J 



-m 1+ ln + --)m 2 



— Jri-m^ — —J rrii 
N 









A7 ml + 


n + - 










A7 ml + 


n+ - 


A?)™ 2 . 



— Jn_m_ — — Jm 2 



N ]m - 



Shi = SB + yn+ — — j JSm 2 — n-JSrri- 
Sh- = ——JSmi — (n + — — J JSm 2 + fn_ — — ) JSm 



sors (2.4) of spins in the two sublattices; because these 
are evaluated for the unperturbed solution (with site- 
independent fields B), there is no need to distinguish 
between \i an d Xi- We thus obtain 



These equations can solved explicitly for the changes 
in the magnetizations; keeping only the leading order 
terms for N — > oo, one finds 



Smi = x+5B 
Sm 2 = -j^x+J [1 



Srri- 



n+X+J - n-X-J\ X+SB 
JjX-J [l - n+x+J - n-X-J] ~* X+SB. 



The 3x3 tensors multiplying SB on the r.h.s. give the 
entries of the susceptibility matrix drrii/dBj (for the 
case j = 1 e I + considered here; as stated above, the 
case j G 7_ can be treated in exactly the same fash- 
ion). We read off that these entries are of 0(1) only 
for i = j, while all off-diagonal terms are 0(1/N). The 
fluctuation-dissipation theorem (AT) then implies that 
all correlations (A<TjA<Tj) between different spins i ^ j 
are 0(1/N); in the thermodynamic limit, mean field the- 
ory therefore becomes exact. 

Appendix B: Direct solution for a special case 

It is possible in some cases to derive the variational 
free energy (2.6) directly from the Hamiltonian (2.1). We 
consider the example of the Heisenberg model (A = 1) 
here to confirm the variational calculations. 

Setting H — 1, the Hamiltonian (2.1) can be written in 
terms of total spin operators as 



H 



4 

'N 



S 2 + + S 2 _ 



1 



2BS Z , 



where we have ignored a trivial constant term. The spin 



operators are defined by 

S± = ^£>*' S = S+ + S-. 

iei± 

Because the quartet {S+, S_, S 2 7 S z } form a set of 
mutually commuting operators, the Hilbert space of 
the iV-particle system is spanned by their simultaneous 
eigenspaces |s+,s_;s,ra). The free energy per spin fx 
is therefore 

•^ = -|? ln E d^- e¥[S+(S++1)+M " +1)1 



2mf3B 



(BT) 



-\s+-s-\ 



Here the sums over s± run from (or 1/2 if N± is odd) to 
N±/2 (= n±N/2). The symbol d 7 ^ denotes the number 
of multiplets of total spin n in a system of M spin-i 
particles. We show in Appendix C that d 1 ^ is given 
explicitly by 



rl n — 



M - 1 
±M + n- 1 



M - 1 
\M + n + 1 



(2n+ 1)M! 



(B-2) 



(\M-n)\{\M + n+\)\ ' 

To evaluate the free energy (B-l), we first carry out 
the sum over m, 



E 



e (2*fl)/3B _ e -(2*fl)/3S 



pB 



-0B 



Setting L — ±(3B, we therefore need to evaluate the 
following sum over s: 

^ e ~^^+l)+{2s+\)L 
s=\s+-s-\ 



= e 



dz _ i »2 
e 2 



2tt 



s+ + s 

E - 

s=|s+-s_ 



2s(L-0/N-izJp/N) 
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X + l 5mi — (^n + — -j-^J J5m,2 + n-JSrri- 



-N J6mi 



n+ -N )J 



^JSm 1 + (n+ - ) J<)rn, 



Sm 2 + n-JSm- 
n- — ) J <Srn 



= SB 

= 
= 0. 



dz 



2ir 



e 2 x 



with 



e 2\s+-s-\{L-P/N-izy/P/N) _ e 2(8 + +s^ + l)(L-f!/N-izy^/N) 
I _ e 2(L-p/N-iz^/JjN) 

Rescaling the integration variable by \/(3N, we thus 
have 



= |s+-s_ 



e (2»fl)/3S _ e -(2*fl)/3S 



e 2/3|a + -s_|(B-l/JV-i«) _ e 2/3(s + + S - + l)(-B-l/iV-iz) 



e -/3B _ e /3(B-2/N-2iz) 
e 20\s+-s-\(-B-l/N-iz) _ e 2f3(s++s- + l)(-B-l/N~iz) 
e /3S _ e /3(-B-2/N-2iz) 



} 



For N — > oo, the third and fourth terms in the curly 
brackets become exponentially small compared to the 
first and second terms, respectively, and can therefore 
be discarded (remember that we assume B > 0). If we 
also make the change of variable z — > iz (which implies 
that the free energy is to be maximized with respect to 
z), we have 



Because n + m + + n-rri- > \n + m + — n_m_|, the first 
term inside the curly brackets becomes negligible for z > 
—B; conversely, for z < —B, the second term can be 
discarded. We are then left with 

1 9 9 9 



The thermodynamic limit of the free energy is then ex- 
pressed as 

f = - lim — In [dz e -^ Nz2 

x E d s N + + d s N _e^ s +( s ++V+ s -( s - +1 V 

f e 2[3\s + -S-\(B-iz) _ e 2f3(s + +S-+l)(B-iz) 

X i e~& B - ej i{ ~ B - 2lz ) 

e 20\s+-s-\(-B-iz) _ e 2f3(s + +s- + l)(-B-iz) 
e /9B _ e l3(-B-2iz) 

The sum over s± can be replaced by an integral in the 
thermodynamic limit: 

s± = -n±m±N, ^2 ^ o n±N / dm± ' 

s± ' ^° 

Also, for large N the expression (B-2) for the combina- 
torial terms d 7 ^ can be simplified to 

1, , s± _1 f N±l 

N ^ d N ± N ^ | [N±{ , _ m±MN±{ i +m±)] , >< 

27V±to± + 1 
JV±(i+m±) + l. 
= n±s(m±) + 0(l/N). 

The free energy thus becomes 

rp i>00 i>\ i> 1 

f = - lim — In / / dm+ / e -0Nf(z,m+,m-) 



f(z, m+,m_) = - n+m+ 

-Tn + s(m + ) — Tn_s(m_) 



2 2 

n_m_ 



■(z + B) x 



|n + m + 



(z > -B) 
(z < -B) 



(B-3) 



Taking the derivative of / with respect to z subsequently 
gives the equations z = — (n + m + + ri-m-) if z > — B 
and z = — \n + m + — n_ra_| if z < — £?. In addition we 
have a local maximum at z = — B if the z-derivative of / 
is negative for z > — _B and positive for z < — _B, which 
translates into the extra solution z = —B appearing for 
|n + m + — ra_m_| < B < n + m + + n_m_. In combination 
the full picture now becomes: 



Elimination of the variable z using the above result 
leads us to a reduced free energy minimization problem 
involving m± only, with 



It is straightforward to verify from Eqs. (B-3), (B-4) 
and (B-5) that the free energy in each phase agrees with 
that given in §3.2. 

Appendix C: Degeneracy factor d 1 ^ 

We derive here the expression (B-2) for the number 
of multiplets of total spin n in a system of M spin-i 
particles. In other words, d 7 ^ is the degeneracy of the 
state \n,n z ), where n is the total spin quantum number 
and n z is any of the allowed values (— n, —n + 1, . . . , n) 
of the z-componcnt of the total spin. 

We proceed by induction over M. For M = 1, we have 
d%j = 1 for n = 1/2 and d%j = otherwise. Now con- 
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f(z, m + , m_) = — z 2 — n 2 + m 2 + — n 2 _m 2 _ — Tn + s(m + ) — Tn_s(m_) 
lim — In ■ 



n^oo N \ smh\3(iz-B)] smh\3(iz-B)] 



gj3N\n+m+ — ri-m- |(— iz— B) g0N(n+m+-\-n—m—)(—iz—B) 

+ ' 



smh[f3(iz+B)] smh[j3(iz+B)} 



f(z, m + , m_) = — ^z 2 — n^m"^ — n 2 _m 2 _ — Tn + s(m + ) — Tn_s(m_ 

^ f gif|n-|_m_i_— TO_m_|(z-t-B) gK"(n + m++n_m_)(2+B) *| 

~ k^oo K ln ( sinh[-/?(.?+B)] + sinh^ + B)] J ' 



sider a multiplet of spin n in an M-spin system. When 
one spin-1/2 particle is added, this multiplet splits into 
exactly one n + and one n — \ multiplet. The excep- 
tion is n — 0, where only a single n = \ multiplet is 
generated. We therefore have 



a M+l 



d 



A I 



+ d 



,n+i 
M 



with the boundary condition 

d n M =0 (n< 



(C-l) 



(C-2) 



Now consider an unbiased random walk with sites num- 
bered by the integers I — 2n + 1 and discrete time 
t = M — 1. The recursion (CT) then tells us that 



-(M- 



a M 



Pt(l), 



(C-3) 



where p t (l) is the site occupation probability at time t 
of a random walk starting from initial position Iq = 2 at 
t = 0. The boundary condition (C-2) simply corresponds 
to an absorbing wall at I = 0. Without this absorbing 
wall, one would have 



Pt(l) = 2" 



in the presence of the wall, one simply subtracts the re- 
flected solution in the usual way to get 



Pt(l) 



L\2(*-*o+*) 



t 



Hi + h + t) 



Inserting this into (C-3) immediately gives the desired 
result (B-2). In writing the above formulae, we use the 
convention that the binomial coefficient (?) is zero when- 
ever k is non-integer or outside the range . . . n. 

Appendix D: Criteria for spinodals, critical and 
tricritical points 

In this appendix, we set out the general criteria that 
we use to find spinodals, critical points and tricritical 
points. Rather than the traditional determinant condi- 
tions due to Gibbs, 7 ' we use a formulation due to Bran- 
nock 8 ) which is more convenient, especially for tricritical 
points. 

Let us assume the free energy per spin is given by 



f(ip), where tp — (ipi...ip n ) is a vector of n (non- 
conserved) order parameters. Thermodynamic phases 
correspond to (local) minima of f(if>) and therefore obey 
the stationary condition Vf(ip) =0. At a spinodal 
point, there is in addition an instability direction Sip 
along which the free energy has zero curvature, implying 
that the gradient of / remains zero to first order: 

(Sip ■ V)Vf(i/>) = (spinodal). (D-l) 

The condition for such a 8xp to exist is 



|M|=0, M = VV/(V). 



(D-2) 



At a critical point, the separation between two neigh- 
bouring (in the ?/>-space) stable phases becomes infinites- 
imal. These phases are separated by an unstable phase (a 
saddle point of /). Constructing a curve ip(s) through 
these three phases (with ip(0) = ip, the state we are 
interested in), we see that (the vector- valued function) 
Vf(ip(s)) vanishes at three infinitesimally separated val- 
ues of s. At the critical point, located at s = 0, these 
three zeros of Vf(ip(s)) coincide, so 

V/(V(s)) = 0(s 3 ). 

Similarly, at a tricritical point we have three stable 
phases coming together, with two unstable phases be- 
tween them, so the corresponding criterion is 

V/(V(*)) - 0(s 5 ). 

Noting that the spinodal criterion (D-l) can be written 
as V/(V>(s)) = C(s 2 ) for the curve ip(s) = ip + sSip, we 
can summarize all three criteria as follows: If there exists 
a curve tp(s) through the point tp (with ip(0) = ip) such 
that 



V/W>(.s)) = 0{s l ) 



(D-3) 



then for I — 2, 3, 5 respectively ip is a spinodal, critical 
and tricritical point. 

To evaluate the criterion (D-3) in practice, we write it 

as 



ds 



V/(V(s)) 







for k = 1. (D-4) 
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B < \n + m + — n_m_ \ : 

\n + m + — n_m_\ < B < n + m + + n_m_ 

B > n + m + + n-m- : 



z = — \n + m + — n-77i-\ 

z = -B . (B-4) 

z = — (n + m + + n_m_) 



(a) Phase A_ : B < \n + m + — n_m_| 

(b) Phase R : \n + m + — n_m_| < B < n + m + + n-m- (B-5) 

(c) Phase A + : B > n + m + + n_m_ 



For the spinodal criterion, this reduces to (D-l) if we 
identify Sip and i/>'(0). For critical points (I = 3) one 
obtains the additional equation 

V((ty • V) 2 /^) + V(V"(0) • V)/(V) = o. 

Taking the scalar product with Sip and using (D-l) elim- 
inates the second term, showing that the criterion for 
critical points is 

(Sip ■ Vff(iP) = (D-5) 

together with (D-l). Following a similar procedure, one 
finds that tricritical points obey the additional condition 

(Sip ■ V) 4 /(V) - 3v • M~ x v = (D-6) 

where 

v = V(5ip-Vff(iP). 

Even though the matrix M has a zero eigenvalue 
(MSip = from (D-l)), the inverse of M in (D-6) is 
well defined: From (D-5), we have Sip ■ v = 0, so that 
v is orthogonal to the corresponding eigenspace. Note 
that while the first term on the l.h.s. of (D-6) is what one 
might have expected naively, the second term cannot be 
neglected: It accounts for the fact that the curve ip(s) 
passing through the three (infinitesimally separated) sta- 
ble phases is generally curved rather than straight. 

We finally note that (D-6) is derived from (D-4) for 
k = 3. In principle, the equation for k = 4 gives an 
additional condition that tricritical points must obey. 
Because of symmetries present in our problem, however, 
this condition is always satisfied in the cases we consider, 
and so we do not give its explicit form. 
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